function y = initial_function1(x)
y = 0.5*(pi-x);
    